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There is no good reason to think that BSM-CP-odd phases(s) will necessarily cause large devi- 
ations in B-physics from predictions of the SM. Therefore, residual theory error in extraction 
of the unitarity triangle can undermine experimental efforts to search for BSM phase(s). We 
stress that final states containing D° or D° in decays of charged and neutral B's can yield all 
the angles of the unitarity triangle with negligble theory error (i.e. O(0.1%)). 

I Introduction & Motivation: why must we target unitarity triangle with "zero" 
theory error 

The two i?-factories have made considerable progress in determining BJoy measurement of the 
^ ■ time dependent CP asymmetry in B°(B°) — > tpK s and related modes™. These machines have 

also been performing remarkably well. Very soon each of them should have ~ 10 8 B-B pairs 
and with improvements in luminosity to > 10 34 cm _2 s~ 1 that are anticipated, 10 9 B-B pairs 
should become accessible in the next few years. In addition, Tevatron experiments CDF/D0, 
and hopefully in the not too distant future -BTeV and LHCB should also enable much larger 
data samples. Furthermore, encouraged by the success of the two B-factories, there is also 
considerable interest in very high luminosities (~ 10 36 cm _2 sec _1 ) facilities, Super-BELLE and 
Super-BABAR. With these anticipated experimental developments it is not enough that we 
can determine (3 with essentially negligible theory error (actually < 1%) we must target a and 
7 extraction also with zero theory errotL Indeed extremely accurate determinations of all 3 
unitarity angles is not just desirable but may well be an essential prerequisite for a successful 
search of the effect of any CP-odd phase in B-decays due to physics beyond the Standard Model. 

In this regard it is important to realize that although the CP-odd phase in the CKM picture 
is of O(l) and not small and it leads to large asymmetries in S-physics, it yields very small CP- 
asymmetries in K decays; recall ex ~ 10~ 3 and e' ~ 10~ 6 even though acp(B — ► ipK s ) ~ 75%. 
There are two important repercussions of this realization. 



1. Failure of the (b — > d) unitarity triangle [UT] due to effects of new physics may well be 
small and subtle. Therefore, residual theory error in the determination of the angles of the 
UT may mask the effect of new physics and thwart experimental attempts to find them. 





Jable 


L: Comparison of s^ 


lilt 


fits. 


Input Quantity 


Atwood & SonE 


Ciuchini et afc 




Hocker et ai 


E> — IT/ , IV , 1 
-fHic = I 'lib/ v cb\ 


nor. 4. ni 7 


.089 ± .009 




n»7 -i- nnfi 4- m a 

.uo 1 in .uuu m .\j±^t 


FB d \jB Bd MeV 


230 ± 50 


230 ± 25 ± 20 




230 ± 28 ± 28 




1.16 ±.08 


1.14 ±.04 ±.05 


1.16 ± .03 ± .05 


£^ 


.86 ±0.15 


.87 ±0.06 ±0.13 


.87 ±.06 ±.13 


Output Quantity 








sin 2/3 


.70 ± .10 


.695 ± .065 




.68 ±.18 


sin 2a 


-.50 ±.32 


-.425 ± .220 






7 


46.2° ±9.1° 


54.85 ±6.0 




56 ± 19 


77 


.30 ± .05 


.316 ± .040 




.34 ±.12 




.25 ± .07 


.22 ± .038 




.22 ±.14 




.185 ± .015 






.19 ±.04 


Am Bs (ps _1 ) 


19.8 ±3.5 


1 7 0+1.5 




24.6 ±9.1 




(2.55 ± .35) x 10" 5 






(2.8 ±.8) x 10" 5 


BR(K + -> tt+vP) 


(0.67 ±0.10) x 10~ 10 






(.74 ±.23) x 10~ 10 


BR(K L -> tt°v&) 


(0.225 ± 0.065) x 10~ 10 






(.27 ±.14) x 10" 10 



2. Search for effect of any BSM phase(s) may well require very large data samples. A model 
independent estimate is very difficult to make. If the asymmetry due to new physics is of 
order ex ~ 10 -3 , then even with a branching ratio of ~ 10 -3 (there are very few relevant 
-B-decay modes that have branching ratios this big) we may need ~ 10 10 B's to find such 
an effect. Efforts at developing the capabilities for large data samples of B 7 s in clean 
environments and/or specialized -B-detectors are therefore very worthwhile. 

These considerations lead us to suggest that while attempts at using proposed methods for a 
and 7 via irn, pn, Kir etc. □ should continue, we should realize that the presence of penguin, 
especially electroweak penguins (EWP), contributions along with the use of flavor symmetries 
could easily lead to theory errors of 10% or more due to any model dependence and theoretical 
assumptions that need he invoked. Clearly on both the theory as well as the experimental front 
methods for determining a, /3, 7 with zero theory error should be vigorously pursued. 

Recall that the current measurements of sin 2/3 ae;ree very well with theoretical expectations; 
the experimental world average sin 2#?£i= 0.78±.offll is completely consistent with theoretical 
expectations, (sin2/?) 5M w 0.70 ± .KB'Bfl. See Table | for a comparison of some theoretical fits. 
It is important though to realize that this test of the SM has serious limitations. A lot of 
theoretical input is used in making these fits as they rely heavily on theoretical evaluations 
of several of the hadronic matrix elements. Improving the accuracy in these calculations is 
extremely difficult and painfully slow. To underscore this we mention two problems with the 
theoretical input being used in these fits. 

First, even for the highly matured calculation of -BxJihere are some reasons to believe that 
the JLQCD (quenched) result, ^^* a s , s ered = 0.86O ± .058H, which has been widely used in the 
past many years, may well be 10 — 15% higher than the true value. This expectation is based on 
results obtained by using the newer discretization, domain wall quarksfDWQ) which has much 
better chiral-flavor symmetry properties. With DWQ both CP-PACS0 and RBC0 get smaller 
values; averaging their numbers one gets Bg WF = 0.758 ± 0.033; (a gain this is quenched). This 
is about ~ 13% below the older JLQCD result. 

Second problem is with the SU(3) breaking ratio, £, which monitors B s versus -Bw oscillations. 
The concern with regard to £ that we have been voicing in the past couple of yearsolli'B is that 
the widely used central value (~ 1.15) provided by some lattice calculations, is quite likely an 
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underestimate and even more importantly the quoted error on ^ of ~ 0.05 appears too low. 
This worry is based on the observation that most lattice calculations of ^ use the "indirect" 
method in which the matrix element of the 4-quark operator is parmetrized in terms of a IB- 
parameter and the pseudoscalar decay constant. In actual numerical calculations of the decay 
constant, the usual practice is to fit linearly to the light quark mass (m q ) dependence. Since the 
orginal matrix element of interest depends quadratically on the decay constant, this procedure 
is unlikely to get the right coefficient^ _fpr example, of m^. Infact, these matrix elements can 
be calculated directly on the lattici3lL3. There is no need to introduce B-paramerters. Use 
of the B-parameters in calculating the AF = 2 mixing matrix elements is purely a historical 
accident. Recall that originally B-parameters were introduced in such calculations for dealing 
with the analogous K — K mixing matrix element. There fx was known experimentally and 
in the phenomenological literature Bk was introduced as a measure of deviation of the matrix 
element for the idealised case of vacuum saturation. For B-mesons, the decay constant is not 
known from experiment and determining that from theory becomes the central issue. One can 
equally well directly calculate (on the lattice) the mixing matrix element without introducing 
decay constant or B-parameteroS. The direct method seems to give largish central value but 
within (rather large) errors is consistent with the indirect method. To be on the safe side one 
ought to use both methods with very good control over errors in each case and then an average 
of the two methods should be used for £. 

Meantime, in light of this observation, we @ adopt a conservative attitude and had used 
£ = L16±0.08 with an error that is considerably bigger compared to Hocker et al. H and Ciuiihini 
et aZ.i; and we have been stressing for quite sometime that their errors are an underestimat£i'0. 
We also examined the effect on sin 2(3 of an increase in £ (including a larger error) along with 
a decrease in Bk, following indications from DWQ. (See Table [2] fromEj). With the size of 
uncertainties currently present, sin 2(3 is hardly effected with {sm2(3) SM = 0.72 ± .10 giving us 
additional confidence in the comparison of the experimental measurements to the predictions of 
the SM. However, the point still remains that use of theory input has its limitations and we 
must try hard to develop methods that can yield angles of the unitarity triangle very "cleanly" , 
i.e. with zero theory error and with no theoretical assumptions. 

We emphasize here that final states containing D°, D° in decays of charged or neutral B's 



can be used very effectively to determine all three angles with essentially zero theory error, i.e. 
<0.1%. 

and B° decays to (e.g.) D°(D°) and K°/K° D°(D°) respectively involve only decays 
via two tree graphs (b — > c and b — > u), no penguin strong or electroweak are involved. Methods 
based on these using direct £3 and time dependents CP respectively can give all three angles 
with no theory error or theory assumptions. In the case of B^ — > K ± D°(D°), the D°(D°) decays 
must proceed to CP-non-eigenstates. Using only D°,D° decays to CP eigenstates as advocated 
in 111 has the difficulty that the suppressed branching ratio is not accessible to experiment as 
D°,D° flavor is difficult to tag in these B-decayf-3. Estimating this branching ratio with the 
use of theory input or assumptions defeats our original goal of determination of the UT with 
zero theory error. A second undesirable feature of using only CP-eigenstates of D° is that the 
resultant CP asymmetry is small ~ (a few %)§; with CPNES methods the asymmetries are 
large. 



1. 1 7 with zero theory error 

This is a uniquely clean method with no theoretical assumption and involving no penguin contri- 
bution, QCD or EW. Interference between two tree graphs, b — > u and b — > c is exploited. Conse- 
quently, the limiting theory error is completely negligible. Furthermore, the interference between 
the amplitudes contributing to common final states of D° , D° that are not CP-eigenstates is large 
resulting in large direct CP asymmetry ~ tens of percents which can be studied, in principle, at 
any B-facilitjc3. Furthermore, while only two modes are essential for the analysis, many modes 
are available. Thus discrete ambiguity in determination of 7 can be removed by use of several 
modes. 

As a specific example one may consider B~ — ► K~D (D ) with D°(D°) — > K + ir~ so the 
overall reaction being studied is just B~ — > K~ K + tc~ . It is important to understand that 
the suppressed branching ratio B~ — > K~D° is not needed and in fact is an output, i.e. is 
determined in the analysis along with 7. The FS (e.g. K~ K + tt~) results from interference 
between two amplitudes one of which is color allowed (B~ — ► K~D°) but doubly-Cabibbo- 
suppressed (D° — ► K + it~) whereas the other is color-suppressed {B~ — > K~D°) but Cabibbo 
allowed {D° -» K+ir~). 

B~ and B + decay amplitudes to two such final states (say, e.g. J2 ± — > iT ± Z) ,Z) with 
D°,D° — > K^ir^ and ET* ± 7r^) result in 4 equations and 4 unknownsli3. The four unknowns 
are the 2 strong phases (one for each final state), the CP-odd weak phase 7 that we are after 
and the branching ratio (denoted by b) of B~ — > K~D° which is extremely difficult to measure 
experimentally due to severe backgrounds. 

One of the advantages of the method is that it can be applied to many modes e.g. B~ — > 
K-,K-*D°,(D°) with D°,D° -> K + tt~, K+p~, K* + tt~, K + a±, K + tt-it + it- etc. Another 
important point is that the method allows to include D°, D° decays to CP-eigenstatesIl! so long 
as one CP-non-eigenstate is also included. So as a specific example one could use 
(i.e. a CPES) with D°,D° -> K + ir~ (a CPNES). The point is that once one CPNES is included 
sufficient number of observables become available to solve for the branching ratio B~ — > K~D° 
as an output along with 70. 

Fig. 1 illustrates use of only two modes, one CPNES (K + ir~) and one CPES (K®ir°) of 
D°,D° assuming A^^ = (# of B-B pairs) x (acceptance) = 10 8 . Solutions to the equations for 
the two modes intersect in four places in the 6-7 plane. Regions with 68%, 90% and 99% CL are 
shown. Multiple solutions are clearly a limitation. Fig. 2 shows the result when several more 
modes of D°,D° are also combined. Now the improvement over Fig. I is significant and 7 with 
a (1-sigma) accuracy of about 7° is obtained; in this calculation true value of 7 is assumed to 
be 60 degrees. 




Figure 1: The likelihood distribution is shown as a function of 7 and b(K*) (which is branching ratio of 
B~ — > K*~ D°) assuming that Nb = 10 s and assuming only the K + ty~ and K s ty° modes of D° are measured. 
The outer edge of the shaded regions correspond to 90% confidence while the inner edge corresponds to 68% 
confidence. The solid lines show the locus of points which give the K + ir~ results while the short dashed curve 

shows the points which give the K b -k° results. 

Figure 2: The likelihood distribution as in Fig.l but now using 6 decay modes of D°,D°. The solution for 
the K + ir~ data is shown with the solid curve; that for the Ksn° is shown with the short dashed curve; K + p~ 
is shown with the long dashed curve; K + a^ is shown with the dash-dot curve; K B p° data is shown with the 
dash-dot-dot curve and the solution for the K* + ir~ data is shown with the dash-dash-dot curve. 



1.2 Extracting j3-a with zero theory error 

The angle 5 = (3 — a + ir = 2/3 + 7, can be obtained by time dependent QP asymmetry 
measurements in B°,B° — » K°D°(D°). As in the case of decays to K^D^^D ) only tree 
graph decays are involved and no penguin, strong or electroweak, enters. Also once again it 
is simply a matter of writing down a bunch of equations and providing enough experimental 
information via observables to render the system completely soluble yielding the unitarity angles 
that we are after. Indeed there is so much redundancy that not only a but also j3 can be 
determined in this way providing a valuable comparison against the (5 determined from B — > 

Some aspects of this method have been previously studied by quite a few authors EHll§!. In 
principle, time-dependent CP asymmetry (TDCPA) measurements in B°(B°) — ► K S D°, K S D° is 
all that is needed to extract 5. However, as in the case of B ± -> K^D , (D°), the L>°, D° flavor 
tagging is problematic. The resolution of this problem using CP-non-eigenstates of D°, D° (just 
as in the case of B^ decays) has been suggested to give information on 5 and also possibly on 
/?§. 

Actually in this case D°(D°) decays to CP eigenstates can also be used. However, if only 
an exclusive CPES (e.g. K s ir°) is used then the number of observables is 3 and the number of 
unknowns, including S, is 4 so not information is available for a separate solution. 

Combining the CPES and CPNES methods seems very effective as the number of unknowns 
involved for the CPES case are just a subset of those for the CPNES case. If one exclusive CPNES 
mode (such as K + ir~) and a CPES mode (e.g. K s ir°) are used then we get 9 observables for 5 
unknowns; if j3 is treated as an unknown even then the system of equations is solvable. 

Of course several CPNES modes can be included. For each CPNES that is included we 
have 6 new observables at the expense of only one additional unknown. Many exclusive modes 
are available e.g. K~ir + (Br ~ 3.8%), K~ p + (10.8%), K*-tt + (5.0%), K*°ir° (3.1%), K*~ p + 
(6.1%), and K*~af (7.1%), for a total of 36%. A nice 4-body mode with all charged track is 




30 60 90 120 150 180 30 60 90 120 150 180 

delta beta 



Figure 3: The Xmin vs - & f° r the toy model calculation given TVs = 10 9 . The thin solid line is the result for 
D° — > K~tt + alone. The dashed line the result for CPES containing Ks together with related CPES containing 
Kl- The dotted lines the result obtained combining K~tt + with CPES containing Ks- The dashed-dotted line 
gives the result for K~ + X alone and the thick solid line combines K~ + X with CPES containing Ks ■ (Note 

S = 110° is assumed here) 

Figure 4: The Xmin vs - P f° r the toy model calculation given Nb = 10 9 using K~ + X with CPES containing 

Ks- (Note P — 25° is assumed.) 



K-TT+TT+TT- (BR > .6%). 

In fact a very nice way to solve for 5 (and (3 simultaneously) is to generalize the above 
exclusive (CPNES) case to inclusive CPNES via D° — > K~ + X. Then the Br~ 53%; one has 6 
observables and 6 unknowns. So it is a solvable system but with discrete ambiguities. A very 
promising way to overcome the ambiguities is to combine this inclusive CPNES case with a 
CPES mode. 

Fig. 3 illustrates the results of our study of extracting 5 via this method. Combining the 
inclusive CPNES (D° — > K~+X) with exclusive CPES seems to do a very good job of eliminating 
the ambiguities and give^L with an error of ±2.5° (the true value of 5 in this case study is 
110°). (See also Table In this example we have used Nb = (number of B-B mesons) 

(acceptance) = 10 9 . 

A very important feature of this method is that you can also use it to solve for (3 very cleanly, 
i.e. no theory assumptions are involved. Infact this method is even cleaner than the B — > tpK® 
method as the latter does receive some (although very small) penguin contributions whereas 
the former has none. However, it is not as efficient as the tpK® method, i.e. more number of 
B mesons are needed to get similar quantitaive accuracy on (5. Fig. 4 illustrates how well the 
method works for determining j3. Here the input used is the one given in the 5th row of Table |3|. 
With Nb = 10 9 the one sigma error on (3 is around 2 degrees. 

Table |] shows a brief summary for determination of the three angles of the unitarity triangle. 
It also highlights the limiting theory error of each method. With the large data samples that 
should become available, it is reasonable to expect that with these B — > K(K*)D° methods all 
three angles could be determined precisely providing an extremely important test of the CKM 
paradigm. A notable feature of these methods is that in charged or neutral B-decays final states 
relevant to extracting all three angles cleanly all contain D° or D°; this should help in improving 
the experimental efficiency. 

Recall that two prominent methods for a, both using time dependent CP, have been studied 
for quite some time. The first methods, requires CP asymmetry measurements in B°(B°) — > 
7r + 7r~ as well as BR for B°,B° — > 7r 7r° and B — > vr^vr . Although these measurements should 
enable one to perform an isospin analysis and remove the penguin contribution the value of a 
thus deduced suffers from contamination from electroweak penguin contribution as EWP evade 
the isospin analysis. To that extent this method for a determination has some residual theory 



Table 3: Attainable one sigma accuracy with various data sets given Nb = 10 9 ; note the 2nd and 5th cases are 

omitted from Fig 3 for clarity. 
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CPNES K-TT+ with K s and with K L 


±5° 
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K~ + X together with K s CPES 


±2.5° 


K~ +X together with Kg as well as 
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error and has to invoke model dependent estimates in the evaluation of the EWP contribution. 

The key experimental difficulty in this set of measurements is the 2tt° mode due to the 
small branching ratio (~ 2 x 10 -7 ) that is expected, made harder by the relative low detection 
efficiency, perhaps also compounded by the fact that the 7r°'s are very energetic giving rise to 
a small opening angle between the 7-pairs. Presumably, these experimental difficulties will be 
surmounted as luminosities improve. Already the two groups have made attempts to measure 
TDCPA in B° — > -k + -k~ which should become quite accurate relatively shortly. However, the 
interpretation of this in terms of the angle a requires important input from theory. 

Another important method wherein isospin analysis, can be used is B — > pit @. In this 
approach one can make use of resonance effects in the Dalitz plot; however, some model depen- 
dence is likely to occur in handling the continuum of B — > 3ir. Once again, EWP are assumed to 
be negligible and this is also an important source of the li mi ting theory error in the pir analysis. 

Recently we have also proposed two other methods for extracting a and 7 that use 
penguin and tree interference effects and therefore are also not completely clean. Table |E| shows 
a sample of these methods that uses penguin-tree interference along with the very clean methods 
of Table ^. It is also important to emphasize that B — > KD° methods require roughly the same 
number of B's as the pir or tttt methods. 

While our emphasis here has been on B + , B~ , B°, B° mesons, if time dependent oscillation 
studies in B S ,B S become feasible, then a clean way to get 7 may also be accessible via final 
states of D£ (D~)K~ (iT + )E2l, or their vector counterparts^. Experimental feasibilty especially 
in a hadronic environment of some of these is studied inE3. 

Let us briefly remark that there is also- a very clean way to get the magnitude of the CP-odd 
CKM phase directly from Kl — > 7t°i^e3'E3. 
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